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SUMMARY 


The eg^uations of motion for the longitudinal dynamics of a tilting 
prop/rotor aircraft are developed. The analysis represents an extension 
of the equations of motion developed in NASA TM X-62,369 to include the 
effects of the longitudinal degrees-of -freedom of the body (pitch, 
heave and horizontal velocity). The development and notation follovr 
that of NASA TM X- 62,369 such that, the effects of body freedom can 
be added to the equations of motion for the flexible vring-propeller 


combination. 
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INTRODUCTION 


This report develops the equation of motion for the longitudinal 
dynamics of a tilting prop rotor aircraft. The effects of v/ing and 
prop rotor flexibility are included in the analysis as well as three 
longitudinal body degrees -of-freedom. The notation and development 
generally follows that of NASA TM X- 62, 369 and is formulated in such 
a way that it modifies that study to include the free longitudinal 
motion of the complete airframe, such that, the influence of the 
wing and prop rotor flexibility on the vehicle dynamics, as related 
to stability and control, can be examined. 

The resulting equations of motion developed in NASA TM X-62,369 
are expressed in matrix notation as follows : 

■^2 -^1 "** -^O^R "^ Ag Cy + AiQ? K Q ~ R B gS 

]? - Cg Xr + G 1 Xr + C q x„ + Cg cf + C-i^ d + C 0 O' + Dgg 

01 - Cx 

. w 

a p X + a n X + a o X = bv + b fvg + Q' F 


( 1 ) 

( 2 ) 

(3) 

(^) 


where x^ is the state vector of the rotor degrees-of-freeomj x^ is the 
state vector of the wing degrees -of-freedom; Vj^ is the input control 
vector; and g, is the gust input vector. P is a column matrix repre- 
senting the forces and moments produced by the prop rotor and a is a 
column matrix giving the linear and angiilar displacement of the prop 
rotor hub due to wing tip motion. 

The development here is concerned with adding the effects of 


vehicle body motion to these equations of motion. The follovring modi- 
fications are developed. Equations ( 3 ) must be modified to include the 
influence of body motion on rotor hub motion. Equations (l) and ( 2 ) 
do not require modification owing to the way in which the equations 
have been formulated. Equation ( 4 ), the wing equations of motion must 
be modified to account for the influence of fuselage motion and in 
addition the fuselage motion equations must be developed. This report 
therefore, is concerned with two items: development of wing/body 

equations of motion and incorporation of the body motion into the hub 
displacement expressions.. The analysis of NASA TM: X-62,369 is not 
discussed in detail in this report. 



MALYTICAL DEVELOBCENT 



A small disturbance approach is used in developing the equations 
of motion. The cehtei’ of gravity of the fuselage is taken as the 
origin of the axis system. X and Z, in general, refer to geometric 
distances measured along and perpendicular to the direction of the 
initial flight velocity. The following specific distances are in- 
volved in the formulation: 

Xr, Zr distance from center of gravity of fuselage to 

spar location at wing root. 

^WL 5 distance from center of gravity of fuselage to 

local effective elastic axis of wing section 
in deflected position. 

Xyj, Zy^T distance from center of gravity of fuselage to 

wing tip effective elastic axis location in 
deflected position. 

Xp , Zp distance from center of gravity of fuselage to 

pylon center of gravity including wing deflection. 
These distances are shown in Figure 1. afc, ,, @o Yo . respectively 
the trimmed flight angle of attack, pitch angle, and flight path angle 
and Vq , is the trimmed flight velocity. 

The fuselage center of gravity motion is specified by perturbation 
velocities and initially along and perpendicular to the triimned 
flight velocity and the rotation of the fuselage is denoted by . 

’ 3 - 


Xf, and Zr are geometric distances, cdaracteristic of the aircraft, 
However, they do depend on the trim angle of attack of the aircraft as 
showhi in Figiu-e 1. 

All of the other quantities depend upon the wing bending and 
torsional deflections. The wing orientation with respect to the 
fuselage is specified by three angles, 6wx, the dihedral, 8wg, the 
wing incidence plus body initial angle of attack (Swg = i^ + ofo) where 
i(^ is the wing incidence with respect to the body reference, and 0*0 
is the trimmed flight angle of attack, and Swg the wing svreep angle. 

In order to obtain results which are a modification of MASA IM X- 
62,369 the wing deflection is expressed in terms of an axis system 
aligned with the wing. 

Assuming that the angles dw^, 6wg, 6W3 are small, as is character- 
istic of this aircraft, the relationship among distances in the fuselage 
coordinate system and the wing coordinate system, where an effective 
sweep angle 6W-, and effective dihedral angle 6w’x introduced 




-X- 


' \ 

- ^ 1 


1 + 6W3 + 6W3 

-X- * 

1 

: Xw 

Y j 

> — 

- 6W3 1 - 8Wx 

•X- 

' < 


z ' 

i 


- 8W2 + 6Wx 1 


! ■ 


The minus sign in front of x is required since x,^ is positive aft in 
Reference 1 . 

-X, 

X = - Xw - 6W3 y^ - 6wa Zw ; ^ ^ 

y = - Xw 6W3 + yw - 6wi (6) 

z = - Xw 6x2 + yw Swx + 


and 6 w* are introduced at this point to reflect the fact that 
al, though the wing is swept, the center spar section is not svrej)t as 
described in Reference 2. Sw’^ and SWg account for an effectiAre 
change in elastic axis position owing to the torsional deflection at 
the root. As developed in Reference 2, it is possible to represent 
this influence on the effective dihedral angle and effective sweep 
angle of the elastic axis as 

6wt = 6 w3 (1 - ?Wb) 


6w- 


6 wi (l - ?Wb) 


where |wb is the torsional deflection at the point where the spar is 
bent, i.e., at the wing root. Thus we have 

x = - Xu - yw 6W3 (l - ?Wb) - 6W3 Zu 

y - Xu 6W3 (1 - 5 wb) + yw - Zw (1 - 5 wb) 

z = - Xu 6 W 2 + yw 6 wi (1 - 5wb) + Zu 


(7) 


( 8 ) 


Therefore, the position of the local wing elastic axis with respect to 
the fuselage center of gravity is expressed in terms of wing coordinates 


measured in the wing axis system. • 

+ X 

= Xr - Xu - yw 6^3 (1 - ?Wb) - 6wg Zu 

(9) 

~ 2r + z 

= Zr + Zu - Xu 6wg + yu 6W3. (1 - §Wb) 

The wing tip deflection is expressed in terms of the same quantities 
with the deflection measured at the tip 



^WT “ " ^WT " yWT ^’'’^3 ” S'Wg Zj^Y 

(10) 

2wt = Z„ + Z^T - Xi^T SWs + YwT Swi (1 - ^Wg) 

The pylon center of gravity is located at a fixed distance from 
the wing tip effective elastic axis location, given hy the following 
distances 

AXp = Zpea 

( 11 ) 

^Zp = 0 

where Zpea is the distance the pylon center of gravity is forvrard of 
^he wing tip elastic axis measured along the wing reference system. 
Consequently 

Xp — X^Y + AXp 

( 12 ) 

Zp = Z^.j Y 

This completes the expression- for the linear displacements of various 
locations of significance. The hub position will be developed later. 

Novr the rotation of each point must be considered. The rotation 
of the fuselage center of gravity is 0^ and consequently the rotation 
Of the wing root is also 0^. Now this rotation must be expressed in 
terms of : the wing coordinate system. This can be done by using the 
inverse of the transformation expressed by (5) 

- A0W1, I ” 1 - 6^3 - Swg I 0 I 

A0„l = 64 1 ■ 0^ (13) 

A'l'wi. ) - 6w^ 1 J ( ^ I 

Therefore 

^0 w L ~ 6vr 3 0 ^ 

: ( 14 ) 



> 


where 0^^^ I'-epresents a rolling of the wing section, and represents 
a yawing of the wing section. The effect of can he neglected by 
assigning that the wing is thin. The complete angular motion of the 
wing section includes the torsional deflection, and is equal to 
0 wu ~ (^f "*■ 0 w) 

At the tip of the wing the rotations are expressed back in the fuselage 
reference frame. The t.*-ing tip rotation consists of torsional deflection 
as well as rotation due to the bending slope. Again using transformation 
( 5 ) 

1 6W*3 6Wg / z\ 

- Sw's 1 - 6w^ I e.^jT 

- 6W2 6Wi 1 J I- xV 

where xV j^epresents the rotation of the tip due to the bending 

slope 



0WT 

+ 6Wg 

Xt 


A0J ~ “ SW3 z j + 

0 W T 

+ 8w'^ 

X T 

(16) 

Aijf-i- ~ “ 6W2 2 T + 

6 w’i 

0WT " 




These quantities are the roll angle, pitch angle and yaw angle of the 
tip expressed in the fuselage axes system. To this must be added the 
rotation, of the. body 6^ so that 

0 T ~ 0p ~ ^^3 ^ T 0WT 6W"^ X-j- (17 ) 

This is also the rotation of the pylon center of gravity. . The above 
expressions describe the linear and angular displacements of the various 
points of interest. 

T. 





Kinetic Enerp;y 

We now proceed to develop the equations of motion for tte body- 
wing-pylon combination using the Lagrangian approach to evaluate the 
inertial terms. 


a. ) Fuselage 

The kinetic energy of the fuselage can be expressed as 

“f “ I “f ^ I 8| 


(18) 


where a reference system travelling at the uniform velocity is used, 

b . ) Wing 

A section of the wing has the following kinetic energy 
y Tw 


KEw = P J ^'5.y„ [(x„ + X - e„); 


y T 

+ 2 + X,, e^)^3 +1 J Iw(e^ + ew)® dy 

+ P J (- 0^1 9f ) <iy 


(19) 


where it has been assumed that the wing is thin such that the roll moment 
of. inertia of the wing section is negligible, 

c. ) Pylon 

The pylon kinetic energy is 
KEp =|mp + 


+ (z^ + Zt + Xp + AXp A^t + (AXp A\^\)®] 

where lo'^ and Ip'^ are the pylon pitch and roll moments of inertia with- 
Vy "^x 

out the rotor measured with respect to the pylon center of gravity 


8 


•o- 


t 


1 


1 


i 


axis location. The pylon inertia in roll is neglected. 

0\v’'ing to tlie manner in vliich the equations were f ox-raulated (Flquations 
) Llironj'ih (M) tin; rotor ol.’l'tK'l,;: ;m'. includocl by ii(:t(:r'5iri.u.lnf; tbrj tifT'^ota 
Oi body laotion on rotor hub motion. Thus, to deterraino the eq.ua bions of 
motion for the ving-body combination the kinetic energy of the system is 


liE = + KEh + KEp 

f . 

Now the modified equations of motion can be developed with the geo- 
metrical considerations described above. Prior to developing these equations 
the modified hub motion equations are developed. 

Hub Motion 

The hub position with respect to the wing tip effective elastic axis 
is expressed by AXh, AZp, again measured in the fuselage reference coordinate 
system. This is the hub location with respect to the fuselage center of 
gravity. The total hub displacement with respect to the center of gravity 
of the fuselage is 


X„=X„t + AXh 

^H~^WT'^AZp, 

where 

AXh = [h C: - (h - h^^)] - [-. Z^^ + h S] 6wg 
AZh = [h C - (h - h^^)] 5wg + [- + h S] 

where , . 

C = cos (6p - Awg) 

S = sin (6p - 6wg) 


( 21 ) 


( 22 ) 



and (h - li , ) and represent the displacement of the elastic axis 
from the wing tip location forward and up due to svreep and dihedral. 
Since (h - h ) and Z are propoi-tional to 6w^ and 6W3 these distances 
are approximately given hy 


AXh = [h C - (h - h^^)] - h S 6W3 

ccti 


AZ 


H = [h C] 6W3 + [- Z^^ + h Si 


The rotation of the hub due to the bending and torsional deflection of 
the vring tip is given by E(iuation (I6) plus the rotation of the body 


The hub displacements arising from these rotations in the direction. of 
the body reference system are 

Axh^p) = - Z^ 0^ - AZh A0t 


^^h/r ~ ( 24) 

Ay Hy^i ~ H 4$ f “ r 


And in addition we have the center of gravity displacements, A x^, Az^ 
with respect to the travelling reference system, and the wing tip dis- 


placements measured in the body coordinate system direction AXjg-ry^B, 

AZ(.j t/s • 

The total hub displacements are 

Ax^ = Ax^ + Axj^ye -Z^ 0^ - AZ^ a 9 t 

Ayn = AZh Aij " AXh Ailf-r ' (25) 

Az h = Az^ + Az -j/g'i" Xh 9^ + AXh A9t 

These must now be resolved into the . hub direction, i.e., along and; 
perpendicular to the shaft . 


10 


In orden to remain consistant with the previous development 


and Azh first resolved to the wing dii-ection 


AXhj^ — 


(26) 


AZhw == - 6w2 AXh 

and then resol\-ed through the angle the prop rotor shaft makes with the 
wing (6p - SWs) , 


Xpjs - Axhw C + Azhw S 


■HW 


HW 


y^Hs ~ AZ^ A?t - AXp, Aijiy 


2hs ^ A2h^j C - AXhw S 


( 27 ) 


The rotations also must be expressed in terms of the shaft axes. Again 
first resolving to the wing direction 


A®jj^ = A®t + Ai|fj 6wg 

A7tw ~ -[• ~ A?t ^^3 

and then in the notation. of Reference 1 


O' ~ Aljfyiig C “ Aij^,| S 
X 


( 28 ) 


a = + A0 t 


Qf^ = A$xw ^ A'I^tw ® 


( 29 ) 


Equations (2?) and (29) express the motion of the hub in terms of wing 
motion and body motion. 

The rotation angles of the wing can now be expressed as,, neglecting 


products of tde wing geometric angles 6wi, awg, 6wg, using Eq.uation (l6) 
in Equation (29). 

n, = z't S - x'-r'C + (6wf C + 6w*3 S) 01^^ 

!S. 

Oy. = zV + 5 Wi Xt + e^-r (30) 

C - xf S +(6wt S - 6w* C) e^T 

The displacement of the hub can be expressed as^ using the relationships 
from Reference 2, pp. 113, for the displacement of the effective elastic 
axis from the tip elastic axis. 

= - yr» 5^1 

^ - ^ea " y-w SW3 
and noting a change in axes such that . 

“ ^HS 

“ ^Hs> Th ~ ^HS ■ 

x„ = (Az^ - 6Wg Ax^) G + (- Ax^ - Az^ 6wg) S 

^wT C + x^ j S t [ (X|^ + AWg-Z k) C + (Zh "S'W'g X|^) S] 6^ 

+ [ h + (h - h^^) G - 2^^ S] + h [- 6W3 z'y + 6w/ x/] 

Yh - h S Zt - Ytw <5^1 Zt + h G Xt - y-rw 6W3 xt 

+ 0„T (- h C 8wi - h S 6wg) : (31) 

Zf. = (Ax^ ,t 6wa Az^:) C + :(az^ - 6W3 Ax^) ,S 

+ [(- 2h + Xh 6W3 ) C + (Xh + Zh 6W2) S] .0^ : V 

- Xwx C + Z - S G + (h^^ - h) S] 0;ij T 


12 



Again these agree with the formulation of Reference 1, with additional 
terms for the effect of fuselage displacement and rotation. 

These results can be expressed in matrix notation as follows, 

O' = c Xj^ + d x^ 

where 



where 

! ^WT “ yxw Iwi 

I - ^WT ~ y TW 1^2 

! ® w. t ~ Pw 

I . 

! and 

I . . ■ ' 

i ' ■ ' . , ■ 

j ^WT “ 'nwCi'^Tw) 'Iw 1 

^WT ~ "nw(yTw) *3.^2 

! The elements of the matrix C are given on page 127 of Reference 1. 

! 

Assuming that 6 Wg is a small angle as has been done in the development, 
some of the terms above can be simplified.. The matrix d which expresses 
the contribution of , the fuselage motion to the hub motion is 


V! 


f 


I ? 
I I 


I 

I 


I 


s a. 


y 


O' 


K 

L 


sin 6p 

cos 6p 

0 
0 
0 
0 


cos 6p 

sin 6p 

0 

0 

0 

0 


1 

Xp, cos Op I 
+ Zh sin 6p I 
Xh sin 6p ? 




/ 



Inertial Foi-ces and Moments 

The kinetic energy of the fuselage, pylon and -wing can be expressed 
as 

= -I- (x| + z|) + I 

yv , . 

+ mdy [(i. + x - 9^.)^ + (z^ + z + 9^)^! 

-7^ 

Yt 

+ IJ* I, (e^ + dy 

-7i 

+ 2 ['I mp [(x^ + Xt - Zp e^)® 

+ + Xp + AXp A0 t ] 

■§■ -^Xp (A'lfi/R + s' AXp (AilfT /l 
2 [ i-Ipy (e^ +A9,)^] + I Ipx Ca^t/r )® 

+ 7 Ipx (A\|fx y{ ) ^ 

The influence of the left and right wings have been included. It is 
assumed that the only wing motion is symmetrical such that 

= TIr / 

Pu = Pr 

and ; . ^ 

6Wit^ = - 6W]_f, 

6w2l “ 9‘^sr 
9^3 u “ " ^^3R. 

It can be seen from Equation (l6) that 

AQt /i.. ~ A8x/ R 

■■ ■ 15 ■ , 


and that 


therefore as can be seen from the kinetic energy expression, the effect 
of both wings is simply to double the pylon and wing terms. 

The generaliced coordinates to be used are the motion of the fuselage 
measured an its center of gravity x^., z^, 6^ and the wing deflection modes 
in two directions and qj^g and the wing torsion p^^. These quantities 
are related to wing deflection by 

~ UtI Iwi 

~ "Hw 1^2 
®w ~ ^ w P w 

where the mode shapes and have been normalized, such that, at the 
■tip T = ^Tw T := 1. The relationship between the wing deflection 

coordinates x and z are given by Equation (8) and the relationships 
between pylon angular motion a 6,, and and the wing deflection and 
torsion are given by Equation (l6). Precisely speaking, the distances 
Xj^L, and Xp and z? are functions of the coordinates, that is the 

wing deflection, however, they can be considered constant, equal to their 
or immed flight value in the following development since the variation in 
these terms will, not contribute any linear terms in the final development.. : 
For reference the time derivatives of Equation (8) and (l6) are repeated here. 

X = - Tljg q.(^3 “ 5^2 % Pwi 

/q\ 

^ ~ “ "Hw 4 iv 2 Swg + tIh qwi 

y is assumed to be constant and does not contribute 
to the kinetic energy 


x-r and z-r are found by replacing by yx^ in the above. 

The rotation rates are 

A§t = - 'Hwt 4wi - ?HT Pw + SWg Tlw'x 4 w2 

A9\ = 6w1 TIwt 4wi + ?HT Pw + 'Hht 4w2 
A\Jt = Swg TljJx 4wi + 5 wt Pw ■ 41wt 4w2 

and 

9^ ~ §w Pw 

The terms of the ec[uations of motion are now evaluated by calculating 

d /SKEx 
dt 

The resulting terms in the eq_uations of motion are listed on the 
following pages in matrix form. The following definitions are employed 

('') = quantity normalized by total aircraft mass, M 

* N 

( ) quantity normalized Py 2 ^b 


Therefore, 


m 


m 


■pz 


m^ + 2m p + 2m p 


■fwp 


M 


mp Z p 


m 


wz 


m 


PX 


m 


wx. 


i o ^yw 

' 'M 

mp Xp 
M 

P^TW 

J , ^WL Mi 
M 
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ro];^ 


nip yxw 


(J]iffldy) yiw 




2 D 




^ea 
y? v-i 


M y^w 


Ipy + mp (Xp 2p) 


M yfw 


fi, dy„+ .fm (x„l ^ z.l) ay , 


M y®w 


JmT] Zi^t, ‘^yw 


1 I 

2 b 


JmTl Xi^;, <iyw 

""■n 


2 ^b 


Jlw 


1 I 

2 b 


lW + ^P 






The fuselage X and Z force equations have been divided by the quantity 
jyQ^R where M is the aircraft mass. The pitching moment equation has been 
divided by M yr^i^ and the wing equations have been divided by ^ 1^. 

The resulting equations of motion for the wing, body are of the 
form - 

Dg f == Evw + g . - 

where 




where Dg has been developed above. Note that in terras of the Nomencla- 
ture of Reference 1, 


D: 


^11 


‘•21 


\ *^12 
! ^3 


where was developed in Reference 1. Dg is given on the accompanying 
page. The following section develops the aerodynamic force and moment 
contributions to the matrices Dj_ and Do- 

In addition, the influence of the rotor forces and moments must be 
added to the body equations. The final equations including effects of wing 
fuselage and rotor are 

Ds \ -f Do = Evw + E g + Hi P + Hg x„ 


vrhere 


i h\ 

I 3-h ?. 

i -d 

I 

— i 

aa 

2C 


P = 


I 

I 

I a O' 


and Hi - 


G-a 

■/V.* 

a 


. i 
== j-al 


2C i 
my I 

a O' 


2C 


mx 


acr 


h is developed in Reference 1 and Gj and Gg are developed in a 
following section.:' 


i 


— In* 


u\T 


■0? H-J 

*” 11%'% pp.. 

+ 'XSu)ifj‘ 


“ 0- TfNu' 


^ •'\ /V I ' 

i ■^'^'n\ L0x I 


'TOp -4-Tl%q^ 


- 't>U3^ tl'\|' 

— t^iQ ti) 




4- n. X, 


Sf> _ 

; Xt,.jqx -r^W2. 

1 W^ j 

! Vr^ 

X"^ nr 

'VCco; 

i +mf' (^e-ZlisSV 


tn's^, 


XviS ^954 
4- S? ^ 


-1 3:^: -2^102 

■ Nvn i ^^'' A ^ 

,-.lSl'sa5'T??|5v.,-t.l\^i^^ 

i-XTN'^ ^W,-n’Sv,3^^ 

-asi iV Swj + n' ^ '5'ui, 


+ n' " ii.. 


Su5 

V- 

4-'T\'’^ W-s-^, *^i- 

“ '^^p 

5£ 

"W^p 


-U^u) -VTCip 


W Xpx 


, '(^ 
“■^yOo Su3 




Winp; AGrodyriamic Forces and Moments 





This section develops the wing forces and moments to be added to 
the wing and fuselage equations of motion. The formulation is compli- 
cated by the fact that the wing has a sweep angle 6 W 3 . The geometry 
of the wing is shown in Figure 4. The lift, drag and pitching moment 
on the wing section can be expressed as 


'L — p Vj^j c a 0^ 


d - I p Vf 3 (C,o + CoQ. 3) 


X 


^ea= “ 2 - P (G, 


AW 


mac Ci 


w 


C.l) 


The total lift, drag and pitching moment acting on the wing is 

L = J ^ dy 

D - J d dy 

M = f m dy 

0 ea 

The normal and chordwise forces, n and c are 
n =<C 
c = d. - t a 

The exciting forces and moments for the wing bending and torsion 
equations are 

W 4: J Z T1 dy 

C = J(d - I a) Vi dy 


22 


f 


noting from Figure 4 that 


c = c cos A 


dy = 


cos A 


We can write 


L = |- p Gw aj,, J Vff a dy 


D - I p cwj Vff (Cpo + 


M = |- p cVos A / (C. 


mac c, , 


Cj dy 


W ■= i p cw J V(f O' T| dy 

C - i p CwJ (Cdo “ (Cu - Cpjj) a) T] dy 


M = I p c^ cos A J (C. 


I _ aw 
mac Cm 


C u ) I dy 


'These last; three expressions are nondimensionalized by dividing by 
N - 

y 2 wnich yields 

= 013 a:„ Jvjf § 


Mc[ws = 012 J (Gdo ' (Cl - Cdj^:) a) 

Ftw 


Mp =021 cos A r (C - Cl) § ^ 

J mac cw ^ y^M 


m 

i - w 

nm ' tt Gf a 


Note that and have been nondimens ionalized by the rotor 
radius K, similarly L, D, and M are normalized by]vr2^R and M yjp’^ Q® 
to yield L'X-, D'^', and MX', 

L* = ^^013 a4 S 

M* = lAl A J 

Prom the geometry of the figure, assigning that the inflow angle, , Is 
small, the velocity normal to the wing leading edge is 

~ [ (^0 cos A - (xi^i + z 5^2 ) “ Z 0^i) 

and the angle of attack of the section is 


a = 0M + 


Vj sin 0^ - [ + Z|,| - Xjj Sjyig] ~ X(^L 0^ 


Vm 


where 

A = 6w 3 + x' 

<^W = (^wa + ) 

,9w = (Sw2 ®f) A + Im P 
Yj = (Vo + sin A : , 


Since many of these quantities are perturbation quantities we may express 
the angle of attack as, retaining only linear terms 
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There are in addition, effects of angular rates, 6t and 0^ accounted for 
^^^ith the following increments to C„ and C 
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We can now determine the terms in the six equations of motion by taking 
derivatives of the above expressions. Denoting the derivatives by 

oMq_i 

'^qixp = 

I 

The following expressions result: 
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The chordwise eguation terms are: 
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The torsion equation terms are 
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These are the terms in the wing eq.uations of motion. 

For the body equations of motion, the lift terms will be similar 

to. the normal force terms applied to the wing with the mode shape not 

present in the integral.' Following a similar notation 

C . = 4o 2C, 0 Vn = C 
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The terms in the x force equation do not involve q.uite the same terns as 
the chordwise wing terms since the resolution is along the fuselage axes 

which are initially aligned vdth the wind. 
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The pitching moment about the fuselage center of gravity is 
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The nondimensionalized coefficients are 
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The various integrals involved in these equations are; 
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Tile control terms for the wing body equations of motion at this point 
involve only the effect of flap delfection. 

The wing equation terms were given in Reference 1 as 


5 “ ^18 ®5 
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The terms in the body equations of motion are 
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Thus the control matrix is 


E = 


J 


The input matrices E and E are modified as follows , Since only 

S 

longitudinal motion is being considered ? lateral gust terms are 
not included. 

b 1 
_£ 

b 

g J 

This matrix is given on the accompanying pages. 

Because of the normalization procedure, in the wing equation the 
derivatives appear multiplied by y and in the body equations by 
Since there are two wings all of the wing terms in the body equations 
appear multiplied by 2. This completes the development of the aero- 
dynamic contributions of the matrices and Dq from the wing. The 
wing stiffness terms carry over directly from Reference 1 along with 
the structural damping and the influence of rotor thrust. 

To be added to these matrices are the influence of the horizontal 
tail, fuselage, gravity forces. 

















Rotor Foi^ce and Moment Contributions 

The Contributions of the rotor forces and moments to the body 
eg^uations of motion are calculated. The rotor forces and moments 
are fox-mulated in a shaft axis system and so the deflections of the 
vring tip must be included in the effect of the rotor forces on the 
body. Linear deflections of the tip are assumed to have a negligible 
effect, however, the influence of the rotation of the tip of the wing 
is included. 

The rotation of the rotor shaft in terms of the torsional de- 
flection of the wing tip and the bending slopes of the wing can be 
expressed as follows: 


1 ^ 
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or 

where 

C = cos p - ) 

S = Sin (6pjp - 6wg) 

This relationship is that given by Equation (30) without the influence 
of body attitude, since the forces and moments are being calculated in a 
body axis system. 

forces and moments in the undeflected system (the 'body axis system) 
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ai-e related to forces and moments in the shaft system which is rotated 


by O' ^5 a , and a by the following relationship 
X y z 
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where 

' X \ ■ may be interpreted as 

Si 

'H 1 


Yo ( 

1 Y ( 


S f 



^S / 

T ) 



The forces and moments relevant to the longitudinal dynamics are 
expressed in terms of the rotor forces in the deflected position as 


X' 

z' 


= H - o Y + O' T 
z y 


- - O' H + O Y +-T 

y X 
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a M + M + o Q 
z X y x^ 


where M is measured -about the rotor hub. The forces are now resolved 

y 

in the body axis directions, parallel and perpendicular to the initial 
velocity direction, and the moment is expressed about the fuselage 
center of gravity 

X = z‘ cos 6p - X^ sin 6p 
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Z = Z* sin 6p + cos 6p 


m'= My + Xh Z'- Z^^ X' 

Now perturbation equations are formulated for the variations in X, 

Z and M. For example 

6X = 6Z^ cos 6p - 6X^ sin 6p 

where 

' 6Z' = - 6cv Ho - ^ 6H + 6a Yq + a 6Y + 6T 

y yo X 0 xo 

and 

6X^ = 6H - 6q' Yo - O' 6Y + 6o' Tn + a 6T 
z Zo - y ° Yo 

where ()q indicates an equilibrium flight value. Proceeding with this 
development the force and moment perturbations applied to the vehicle 
by the rotor can be expressed as: 

Pb - G-i F + Gg 


where 






when 
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Pw 

The matrices and are given on the accompanying pages. The following 
notation is used 

C' = cos 6p 
= sin 6p 

and the subscript zero is dropped with the understanding that all matrix 
elements are evaluated at the trim flight condition. Both rotors are 
accounted for in the matrix elements. 




















Gravity Terms 


Also the effect of the weight of the aircraft must he added in 
the fuselage equations. The forces and moments to be added are 

Z = - W cos (Yo + 

X = W sin (yo + Sf ) 

M = W [ZcG sin (yo + 6^ ) " Xcq cos (Yq + 6^) 

wh4re Yo is the initial flight path angle and Xcq > Zcg is the location 
of the total aircraft center of gravity from the fuselage center of 
gravity. Assuming that 6^ is a small angle these forces and moments 
can be written in perturbation form as 

Z - - W cos Yo + W sin Yo 6^ 

X = W sin Yo + W cos Yo 

M ^ W (ZcG sin Yo "Xcg / cos Yo) 

+ W (ZcG cos Yo + Xce sin Yo ) 

Fuselage Aerodynamics 

The fuselage is assimaed to experience only a drag force and a 
pitching moment which can be expressed in normalized form as 

The perturbation terms for the equations of motion would be 
6D* = 2 0og % C 0 f 6 V + 0O2 ; p 

6M* = S <ho Vo 6V Sj 010 Vo= 6a, 

C/ 

k3 ' ' 


and 
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Hoi’izon.tal Tall Aerodynamics 

The horizontal tail is taken to be located a distance -t-r behind 

the fuselage center of gravity and a distance above the fuselage 
> 

center of gravity. , Only the lift of the horizontal tail is considered 
as contributing to the motion eq^uations. The drag of the tail surfaces 
is included with the drag of the fuselage. The contributions to the 
vertical force and the pitching moment are 
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The lift of the horizontal tall suitably normalized may be expressed as 
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The next to the last term accounts for the downwash lag at the tail, 
and T is the elevator effectiveness,. . As a consequence we have the 
folloiving terms in perturbation form. 
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The pitching moment derivatives follow directly from these expressions. 
The pitching moment derivatives are related to the lift derivatives by 


Ki ) 'I 


H T 


Vt W H T 


The trim tail lift coefficient is . 

G(_Tq = a j (ij + ao - g + t ) 

The dcrs'fnwash model used here is satisfactory for forward flight but is 
probably too simple for the low speed condition where the downwash con- 
ditions at the tail are q.uite complex. 

The matrix terms comprising ADxj ADq and AE are given on the 

S 

following pages. (C ) contributes to the matrix Dg. 

ZZ rt 
^ HT 
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EQUATIONS OF MOTION 


The final form of the "wing hody equations of motion are 

Do X + D-f X + Do' X = Ev + E'g + E.6 + H, P + Ho x 

3 V ^ V ° V 'w g 6 e ^ ^ V 

The matrices , Dq and e' include the effects of the fuselage, hori- 

g 

zontal tail, and the qravity forces and moments. E^ is the influence 
of the elevator. Thus 


d'i = Di + ADi 
Uo = Dg + ADg 



The matrices AD 1 , ADg and AE are given: on the accompanying pages. 

o 

The equations relating hub motion & to wing motion x'yj and x^ are 

O' = c X m + d x^ 

These two equation's replace the lower two sets of equations given 
on page l4l of Reference 1 for the dynamic motion of the vehicle 
with the fuselage free to move in the longitudinal plane. 
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The complete set of equations of motion are: 
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TRIM CALCULATIONS 


The ec^uations which determine the eq_uilibrium flight condition of 
the aircraft in normalized form may be written as follows: 
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If the eg^uilibrim deflections of the -wing are included in the 
! 

trim calculation then the solution for trim becomes q.uite complex. It 


seems unlikely that the equilibrum deflections will have a significant 
influence on the trim of the aircraft and therefore they will be neglected 
in the solution for trim. With this assumption, 
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The simplest trim problem is the case in which the rotor shaft is 
aligned v/'ith the initial velocity, that is, when 6p = 0, or in other 
words when ip = ^ (ij^ + cVq ). In this case, it is. further assumed that 
cyclic pitch is not used for control such that also in trim and 9xs = 0 
then the rotor trim condition is a perfectly axisymmetric case and 


52 


and 
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2C 
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= 0 


Given the geometry of the aircraft, and selecting the initial 
flight velocity Vq (actually the advance ratio) and the flight path 
angle Yqj 'blie three equilibrium equations may be solved to determine 
the trim values of the airplane angle of attack o'!, , the elevator angle 

2Cj 

6 , and the rotor thrust, . 

e’ a cr 

For the more general case in which 6p is not equal to zero, but 

/ 

is Still a small angle, that is, when airplane flight is being ‘ 
considered, linearized expressions are employed to calculate the rotor 
inplane force and the rotor pitching moment. 

For this more complex case in which the blades are assumed 
to be torsionally rigid and no cyclic is applied for control, the 
following equations are involved. The gimbal motion is determined from 
(I* (vf - 1) + Kp Y Mpi) Pac = he + 3gs 

(ij {vi - 1) + Kp Y Mpi) gee = he + yy^ hs + Y Sp 


2Ch 

— = Y ^ (H^ + R^) 6p + (y Rp - yH^ - YKp Rpi) he 

+ (YHp+ YSp - Y Kp Hpi) Pgs 


The aerodynamic derivatives in the; above expressions do depend upon the 
trim thrust to some degree. However, it would be expected that the effects 
of rotor inplane force have only a small effect on the force equilibrium. 
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The solution for trim can therefore proceed as in the purely axial 
flow case with the selcted small incidence of the shaft. Once the 
angle of attack of the aircraft is computed such that an initial 
value of 6p is obtained, the flapping coefficients can be calculated 
and consequently the equilibrium values of the inplane force and the 
hub moment can be calculated and the trim calcifLation repeated to 
account for these effects. 

Once this procedure is completed, expressions from Reference 1 
/ - 

can be used to calculate the trim values of the remaining rotor forces 
and moments which are used in the equations of motion. 

If rotor cyclic is introduced into the trim calculation, then the 
relationship between cyclic and elevator angle must be selected. 
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